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Low energy properties of half-filled Fibonacci Hubbard 
models are studied by weak coupling renormalization group 
and density matrix renormalization group method. In the 
case of diagonal modulation, weak Coulomb repulsion is irrel- 
evant and the system behaves as a free Fibonacci chain, while 
for strong Coulomb repulsion, the charge sector is a Mott in- 
sulator and the spin sector behaves as a uniform Heisenberg 
antiferromagnetic chain. The off-diagonal modulation always 
drives the charge sector to a Mott insulator and the spin sec- 
tor to a Fibonacci antiferromagnetic Heisenberg chain. 
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Since the discovery of high Tc oxide superconductors 
and heavy fermion materials, the strongly correlated elec- 
tron system has been the most important subject of re- 
cent condensed matter physics. Even in the insulating 
phase, the quantum magnetism in these and related ma- 
terials have been attracting the wide interest from theory 
and experiment M . Another remarkable finding in recent 
condensed matter physics is the discovery of quasicrys- 
tals iQ. The electronic states in quasicrystals are not 
trivial even in the simplest case of one dimensional free 
fermions. For the Fibonacci lattice, the beautiful multi- 
fractal structure of the single particle spectrum and the 
wave function have been revealed by means of the renor- 
malization group method |^,Q . 

Nevertheless the interplay between the quasiperiodic- 
ity and strong correlation in quantum magnetism has 
been rarely studied except for the recent bosonization 0] 
and density matrix renormalization group (DMRG) stud- 
ies for one-dimensional Heisenberg chains. M Experimen- 
tally, several kinds of quasicrystals with local magnetic 
moments have been synthesized recently M. In this re- 
spect, the quantum magnetism in quasiperiodic systems 
must be a promising field in the condensed matter physic 
of next decade. 

Although one dimensional quasiperiodic magnetic ma- 
terial is not available at present, one of the plausible can- 
didate would be an artificial structure such as quantum 
dot array Iq]. In such realistic situations, the iteneracy of 
electrons becomes important. For the realization of one 
dimensional Fibonacci antiferromagnet, it is necessary to 
specify the parameter regime in which quasiperiodicity 
comes into play. In the present work, we therefore investi- 
gate the Fibonacci Hubbard model in which the coupling 
between spin and charge degrees of freedom produces a 



rich variety of ground states even in the half-filled case. 
Our Hamiltionan is given by. 
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for off-diagonal modulation. The superfices d and o rep- 
resent the diagonal and off-diagonal modulations, respec- 
tively. The operators aj ^ and Oi^o- are creation and ani- 
hilation operators of fermions with spin a{=^ or I) and 
'T-i,CT = CLi^ai^a-- The open boundary condition is as- 
sumed. The on-site coulomb interaction is denoted by U. 
The transfer integral ta- 's (= t^ or ig) or the on-site po- 
tential Va^'s {— Vj^ or Vq) follow the Fibonacci sequence 
generated by the substitution rule A -^ AB, B —^ A. 
The modulation amplitudes are defined by Ai = ^A ^ ^B 
and Ay — V^ — Vg. For the off-diagonal modula- 
tion case, the average transfer integral t is defined by 
t = {tpj) + t]3)/(l + (t>) where (j) is the golden ratio 
= (1 + \/5)/2. In the rest of this paper, we concen- 
trate on the half-filled case. 

First, we employ the bosonization method in the weak 
coupling limit [/, | At |, | /S.V |<< t, to obtain the follow- 
ing bosonized Hamiltonian 
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with 



2/^ = Ua/irvF^ vp 
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Ua — w_F-\/l ~ Ua/nvp, Up = vp^/l+lJa/^^up 

The boson fields (jfp and (/)„ represent the charge and spin 
degrees of freedom, respectively. The momentum densi- 
ties conjugate to them are denoted by lip and IIo-. The 
lattice constant, fermi velocity are denoted by a and vp. 
The ultraviolet cut-off denoted by a is of the order of a. 
The function W{x) represents the Fibonacci modulation 
of amplitude unity. 

Following Vidal et al. Isj , we obtain the weak coupling 
renormalization group (WCRG) equations for the cou- 
pling constants by the standard technique ll^ as. 
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where yq{0) = aXW{q)/vF with A = Ay/2 for the diag- 
onal modulation and A = At for the off-diagonal mod- 
ulation. The Fourier components of M^fx) is denoted 
by W{q) whose explicit form is given in H. The renor- 
malized short distance cut-off is given by a{l) — aeK In 
eq.(|^), the summation over q is performed for q — 2TTm/n 
with m € [1, n — 1] where n is the generation of the 
Fibonacci sequence and R is the Gaussian ultraviolet 
regulator R{x) — e~^ . The -|-(— ) sign in eq.(H) is for 
the off-diagonal (diagonal) modulation case. The cor- 
rections to the velocities Up and Ua are neglected since 
they give the higher order corrections to the renormaliza- 
tion of other quantities. Similar set of equations for the 
alternating potential are derived by Tsuchiizu and Suzu- 
mura fl^ . It should be noted that the scaling dimension 
of the quasiperiodic modulation term is non-trivial be- 
cause of the presence of self-similar function W{x), so 
that we have to resort to numerical calculation to solve 
the WCRG equations. 

We have also carried out the numerical calculation us- 
ing the DMRG method to obtain insight into the strong 
coupling regime which is inaccessible by the WCRG cal- 
culation. In the numerical calculation, we take t = 1 



to fix the energy unit. To reveal the bulk properties 
of the Fibonacci chains, it is useful to investigate the 
behavior of the average of physical quantities over all 
possible finite length subsequences of infinite Fibonacci 
chains as discussed in |^,|l^ . It should be also noted that 
the number of the n-mcmbered subsequence is equal to 
n + 1 |lj|. The number m of the states kept on each 
DMRG step ranged from 120 to 300 depending on the 
values of parameters. The convergence with respect to 
m is checked. If weak m-dependence remains around 
m = 300, the m-extrapolation is carried out for each 
energy eigenvalue E{m) using the extrapolation formula 
E{m) ~ ^(oo) -f- c/m'^ |l||. 
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FIG. 1. The Af-dependence of < InA^.c > for the Fi- 
bonacci Hubbard model with f/ = 4 for diagonal modula- 
tion with AV = 2.0 (squares) and off-diagonal modulation 
with At = 0.4 (circles). Filled (open) symbols represent the 
charge (spin) gap . 

Numerical solutions of the WCRG equations are clas- 
sified into two categories according to the behavior of j/p. 
In the following, we discuss these two regimes separately. 
{a,) Strong coupling regime 

Hyp increases in the course of renormalization, it grows 
more rapidly than other parameters, so that the phase (/)p 
is fixed to 7r/2\/2. The growth of the Fibonacci modula- 
tion term G is suppressed as expected from — j/p term in 
the RHS of eq. (@). In this case, the charge gap opens 
and the ground state is a Mott insulator. For diago- 
nal modulation, the Fibonacci modulation term, which 
is proportional to cos\/2(j)p, vanishes. Thus the low en- 
ergy behavior of the spin sector is described as an anti- 
ferromagnetic uniform Heisenberg chain. On the other 
hand, for off-diagonal modulation, the Fibonacci mod- 
ulation term is proportional to sin v20p which is fixed 
to 1. Therefore the spin sector is renornialized to the 
Fibonacci antiferromagnetic Heisenberg chain whose be- 
havior is discussed in detail in MM. It should be noted 
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that Kp is always less than unity because ^p(O) < 1 and 
< 0. For the off diagonal modulation, therefore, 
-^ is always positive and the above behavior is always 
realized. 



Typical cases are numerically demonstrated by the 
DMRG method in Fig. || with C/ = 4 for the off-diagonal 
modulation with At = 0.4 and diagonal modulation with 
AV = 2.0. The system size dependence of the average 
of the logarithm of the spin gap As and charge gap Ac 
are shown. It is verified that the charge gap tends to a 
finite value as iV — *■ oo for both cases. For the diagonal 
modulation, the spin gap behaves as < In As >~ — IniV 
with slope unity which is the Luttinger liquid behavior of 
the uniform Heisenberg chain. For the off-diagonal mod- 
ulation, the size dependence of the spin gap is well fitted 
by the formula < In As >^ —N'^ as in the Fibonacci 
Heisenberg chain ||] . It should be noted that no trace of 
Fibonacci modulation remains in the size dependence of 
the spin gaps of the diagonal case even though the mod- 
ulation amplitude is 5 times larger than the off-diagonal 
case. 

In the limit of strong U » t^, tg, our Hubbard model 
can be mapped onto the Heisenberg model as. 
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where Si's are the spin 1/2 operators. In the case of off- 
diagonal modulation the exchange couplings Jq/s (= J^ 
or Jg) follow the Fibonacci sequence as, J^ — 2t\/U 
and Jg = 2try/U. On the other hand, for the diago- 
nal modulation, the exchange coupling is determined by 
U and Va of the sites on the both ends of the bond. 
Therefore the exchange coupling can be indexed by the 
pair of letters which appear in Fibonacci sequence as 
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Jaa = ^i /U and Jab = -^BA 

t^ /{U + Vg — V^). The sequence BB does not appear. 
We call this type of modulation as pairwise Fibonacci 
modulation. The system size dependence of the energy 
gap calculated by the DMRG method for both types of 
antiferromagnetic Heisenberg chains are shown in Fig. g 
with A J =1 Jp^^J^ I "-"^ I '^AA~'^AB I- "^^^ energy unit 

is J=(0JA + JB)/(l + '^)or(0JAA + "^AB)/(l + '^)- ^s 
expected, the spin excitations scales as As ^ exp(— cA^"^) 

for the Fibonacci Heisenberg chain and as As ~ — In A^ 
for the pairwise Fibonacci chain, 
(b) Asymptotically free regime 

For the diagonal modulation case, the competition be- 
tween the quasiperiodic modulation and the Coulomb in- 
teraction can take place for small U . In this case, both 
yp and y^ are renormalizcd to negative values and the Fi- 
bonacci modulation term G grows under renormalization, 
so that none of them arc dominating. Typical example 
of the numerical solution of WCRG equations is shown 
in Fig. |(a) for AV/t = 0.4 and U/2iit = 0.02. Com- 
paring this behavior with that of the free fermions with 
Fibonacci potential shown in Fig. H(b) with /SV/t = 0.4, 
we find both flows are quite similar. This implies that 
the main contribution to the renormalizcd quantities are 
generated by the back scattering due to the Fibonacci 



modulation and the bare Coulomb interaction does not 
make essential contribution. Therefore we conclude that 
the weak Coulomb interaction is irrelevant and the low 
energy spectrum scales as < In As,c >^ ~z In A^ similarly 
to the free case. The numerical results by DMRG shown 
in Fig. ^(a) for U = 0.4 and AV" =1.2 also support this 
conclusion if compared with the exact diagonalization re- 
sults for the free case with AV =1.2 shown in the same 
figure. For comparison, Fig. H(b) shows the same quan- 
tities for the off-diagonal modulation with U = 0.4 and 
At = 1.2. In contrast to the diagonal modulation case, 
there is a clear evidence of the finite charge gap. The spin 
gap decreases too rapidly to estimate the precise value for 
large systems A^ > 36. It is, however, consistent with the 
Fibonacci Heisenberg type behavior < In As >~ —N'^ 
H rather than the power law < In As >~ — zlnA^. 
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FIG. 2. The A-dependence of < InA^ > for the Heisen- 
berg model with (a) Fibonacci pairwise modulation and (b) 
Fibonacci modulation with J = 1.0. 
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FIG. 3. The i-dependence of the renormalized parameters 
with AV/i = 0.4 with (a) U/2-wt = 0.02 and (b) f/ = 0. The 
ultraviolet cut-off a is taken equal to a. 
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FIG. 4. The A^-dependence of < In As,c > for the Fibonacci 
Hubbard model with U = 0.4 and f/ = for (a) diagonal 
modulation with AV — 1.2 and (b) off-diagonal modulation 
with At = 1.2. Filled (open) symbols represent the charge 
(spin) gap . 

Following the argument of y , the low temperature be- 
havior of the thermodynamic quantities are deduced from 
the above scaling behavior of the low energy spectrum. 
For the off-diagonal modultion, the magnetic suscepti- 
bility X behaves as x ~ l/(T(lnr)^/'^) and magnetic 
specific heat C as C - l/{\iiT)^+^/'^. This also holds 
for the Fibonacci Heisenberg chains. For the diagonal 
modulation with large U, we expect the Luttinger liquid 
behavior x ~ const, and C ^ T, while for the diagonal 
modulation with small U, we expect the free Fibonacci 
chain behavior x '^ l/T^/''^ and C - l/T^/\ 

In summary, we find that the Fibonacci repulsive Hub- 
bard model at half-filling shows a variety of ground states 
depending on the types and strength of modulation in 
contrast with the free fermion Fibonacci chains which is 
critical irrespective of the type or strength of modulation. 
For the off-diagonal modulation, the effect of Coulomb 
interaction is most drastic. The ground state is always a 
Mott insulator and the spin sector behaves as an antifer- 
romagnetic Fibonacci Heisenberg chain. On the contrary, 
for the diagonal modulation, both spin and charge sectors 
behave as free Fibonacci chains if the Coulomb interac- 
tion is weak enough. This implies that the conventional 
free theories ^ are approximately applicable to the di- 
agonal modulation case as far as the electron-electron 
interaction is weak. This is in contrast to the case of 
spinless fermion chains in which nearest neighbour repul- 
sive Coulomb interaction is always relevant [Hj|] for both 
diagonal and off-idagonal modulations. Even in the diag- 
onal modulation case, the ground state becomes a Mott 
insulator if the Coulomb interaction is strong enough. In 
this case, however, the effective exchange modulation in 
the spin sector is irrelevant and spin sector behaves as an 
antiferromagnetic uniform Heisenberg chain. 

Our calculation suggests how to realize different types 
of Fibonacci electronic system using quantum dot arrays. 
To realize nearly free Fibonacci chain, the local poten- 
tial of dots should be modulated and the charging en- 
ergy should be reduced. On the other hand, to realize 



Fibonacci Heisenberg antiferromagnet, it is the transfer 
integrals (or distances) between the dots which should be 
modulated. 
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